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Abstract. Comparisons of particle production from high-energy ion collisions with next-to-leading order 
perturbative QCD calculations show good agreement down to moderate transverse momentum values. Dis- 
tributions of azimuthal angle differences between coincident hadrons in these collisions support a partonic 
origin to the particle production, again down to moderate transverse momentum values. The rapidity de- 
pendence of inclusive and coincident particle production can therefore be used to probe parton distribution 
functions down to small momentum fractions where theory anticipates that parton saturation could be 
present. This paper describes how such experiments could be completed. 

PACS. 12.38 Qk - I3.88.+e - 24.85. +p 



1 Introduction 

Comprehensive measurements of p+p, d+Au and Au+Au 
interactions at ^/snn— 200 GeV by the RHIC experiments 
strongly suggest that the central collisions of two gold nu- 
clei lead to a new form of matter that appears opaque to 
high transverse momentum (j>t) hadrons pp. This dense 
matter evolves from an initial state produced by the col- 
lisions of the low- a; gluon fields of each nucleus |2J. Un- 
derstanding this initial state is the first topic mentioned 
in a recent report on scientific opportunities in heavy-ion 
physics. "Upgraded forward instrumentation" was iden- 
tified as needed to elucidate the properties of the initial 
state 0. Nucleon gluon density distributions are deter- 
mined by global fits to data @JS], but the low-x nuclear 
gluon distribution is not yet known |fj][7). The nuclear 
gluon field distribution might be naively expected to result 
from a convolution of the gluon density distributions of all 
the individual nucleons. However, there is indirect experi- 
mental evidence from RHIC that the low- a; gluon distri- 
bution in a large nucleus like gold is reduced, or shadowed, 
from the nominal superposition of the distributions of the 
included protons and neutrons, a phenomenon described 
as saturation. In this document we describe a possible 
measurement of the gluon distribution in a large nucleus. 

In the early runs at RHIC, we have demonstrated that 
forward electromagnetic calorimeters (Fig. [TJ, the For- 
ward 7T° Detector (FPD), can be used to measure ir° pro- 
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Fig. 1. Layout of the Solenoidal Tracker at RHIC 



duced at large pseudorapidity (rj = —In tan 9/2, where 
9 is the polar angle of the produced particle) in p+p 
and d+Au collisions at ^7^=200 GeV at STAR [TU] , 
We have recently proposed to assemble a Forward Meson 
Spectrometer (FMS) that will be operated during future 
RHIC running periods by the STAR collaboration to en- 
able measurement of the gluon distribution, xg{x), in nu- 
clei in the range 0.001 < x < 0.1. The function g(x) gives 
the differential probability to find gluons with a fraction 
x of the longitudinal momentum of the parent nucleon 
(Fig. I2J). The FMS will cover the range 2.5 < 77 < 4.0 
and give STAR nearly hermetic electromagnetic coverage 
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Fig. 2. (Left) The gluon distribution in the proton [I]. Note 
the rapid rise in xg(x) for x < 0.01, a discovery made at 
HERA based on studies of deep inelastic scattering (DIS), us- 
ing electron(positron)+proton collisions at -^=300 GeV [111 
I12| . (Right) Values of DIS kinematic variables x and Q 2 where 
nuclear data from fixed-target DIS experiments at much lower 
y/s constrain the nuclear gluon density for x > 0.02 |5I7| . 



in the range — 1 < rj < 4. The FMS will allow correlation 
measurements between forward mesons and photons with 
signals from the full STAR detector, including the barrel 
and endcap electomagnetic calorimeters (BEMC,EEMC) 
and the forward and midrapidity time-projection cham- 
bers (TPC). Exploiting the capabilities of RHIC and the 
existing STAR detector, and assuming simple 2-body kine- 
matics, the FMS will allow measurement of the gluon den- 
sity in protons and in nuclei down to x ~ 0.001. 

With the addition of the FMS, which has > 25 times 
larger arcal coverage than the FPD we view as its proto- 
type, STAR will be able to achieve at least three important 
and new physics objectives: 

— A measurement of the gluon density distributions in 
gold nuclei for 0.001 < x < 0.1, thereby extending 
our current knowledge and including an overlap region 
that tests the universality of the gluon distribution. 

— Characterization of correlated pion cross sections as a 
function of px to search for the onset of gluon satura- 
tion effects associated with macroscopic gluon fields. 

— Measurements with transversely polarized protons that 
are expected to resolve the origin of the large trans- 



verse spin asymmetries in pj 
for forward it production. 
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In d+Au collisions, the FMS will face the deuteron 
beam and will see neutral pions produced by large- a; quarks 
in the deuteron interacting with the low-x gluons in the Au 
nucleus. The key to this analysis is the detection of a sec- 
ond particle in coincidence with a triggering particle in the 
FMS. The coincident signal might be a high-py track or jet 
detected in the TPCs or it might be a 7 or n° detected in 
the BEMC or EEMC. For x < 0.01, the coincident particle 
will be a second 7/77° detected in the FMS, whose large 
acceptance makes this coincidence measurement possible. 
The spatial dependence of the nuclear gluon density [Hfl 
I14j will be investigated by analyzing two-particle correla- 
tions as a function of particle multiplicity in the Au beam 
direction with the existing STAR subsystems. 

Analysis of the kinematics of the relative momentum 
between the trigger particle and the coincident particle 



allows us to determine g(x) in the gold nucleus under 
the simple assumption of elastic scattering of collinear 
initial-state partons. This measurement of the gluon den- 
sity provides the essential input to the simulation codes 
that attempt to determine the energy density achieved 
when heavy nuclei collide, in the state which could ex- 
pand to become the quark gluon plasma. 

The same correlated particle analysis will allow us to 
study the physics of the parton saturation region, if it ex- 
ists for Q 2 < 4 GeV 2 . This physics is associated with the 
transformation from a parton-dominated picture of the 
nuclear gluon distribution to a picture for which macro- 
scopic QCD fields might play a role or provide the most 
approporiate physics description. The FMS granularity 
will enable measurement of the azimuthal angle (4>) of 
the trigger pion and a coincident pion. The peak in the 
A(j) = cjy^i — distribution at 180°, the classic signature 
of parton clastic scattering, is expected to broaden [T5] 
or disappear |lfi| (ie., the forward jet becomes a monojet 
and the recoil hadron rapidity distribution is modified) 
when scattering from macroscopic gluon fields dominates 
scattering from single gluons. 

The analysis of FMS-triggered events at STAR will 
also be used in polarized proton running where the ex- 
tended calorimetry acceptance will greatly enhance our 
ability to determine how quark and gluon fields conspire to 
share the proton's | unit of spin. Polarized deep inelastic 
scattering (DIS) experiments found that the intrinsic spin 
of quarks and antiquarks contribute only ~ 20% to the 
nucleon spin, contrary to early theoretical expectations of 
> 60%. A prime objective of the RHIC spin program is 
to understand how gluon spin and parton orbital angu- 
lar momentum play a role in this "missing spin puzzle" . 
The correlated pion analysis of the FMS and the analy- 
sis channels it opens will play a crucial role in answering 
these questions and could lead to resolution of the long- 
standing question about the origin of the large transverse 
single spin asymmetry in forward pion production. 



2 Nuclear gluon densities 

A central objective in high energy physics has been the 
systematic characterization of parton (quark and gluon) 
density distributions 0fl5] . As a consequence of factoriza- 
tion theorems we know that there is a class of high-p^ 
two-parton (leading twist) experiments that can be under- 
stood in terms of an initial state of independent partons 
within a proton. Within this framework, the part of the 
cross section due to a particular sub-process is equal to 
the product of a calculable parton-level cross section and 
the two universal initial state parton probability densities 

0" \Xn : OC CT nm fji (x n ) fm (Xm ) . 

The parton densities f(x) are universal properties of 
the proton, applicable in all hard scattering processes, and 
in most cases (but not all) refer to the positive definite 
probability density to find a parton "n" (n = q for quark, 
q for antiquark and g for gluon, with f g (x) referred to as 
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g(x)) carrying a fraction a x" of the parent nucleon lon- 
gitudinal momentum; a; is a kinematic variable in DIS. 
We combine contributions from all partonic sub-processes 
that lead to the same final state and account for contribu- 
tions that come from all possible x values by adding the 
partial cross sections. 

The nucleon gluon distribution xg(x) is known in the 
region 0.001 < x < 0.01 (Fig. I2J) but the nuclear distribu- 
tions are not. Our present understanding of how parton 
distribution functions (PDFs) are changed when nucleons 
are bound in a nucleus is primarily derived from DIS of 
charged leptons from nuclear targets. The charged lep- 
tons used in DIS interact with the electrically charged 
quarks, not with the electrically neutral gluons, and pro- 
vide measurements of structure functions, Fi. In the par- 
ton model, F 2 (x,Q 2 ) = x*£ n el[q n (x,Q 2 ) + q n (x,Q 2 )}, 
where is the squared electric charge of the quark of 
type n and Q 2 is the squared four momentum transfer of 
the scattered lepton, equated to the square of the scale at 
which the parton substructures are observed. Gluon den- 
sities are determined from the QCD evolution equations 
[T7) applied to scaling violations of F 2 measured for the 
nucleon over an extremely large x and Q 2 range at the 
HERA collider [TlHr2] . Sensitivity to g(x) in DIS is ap- 
proximately given by the Q 2 variation of F 2 at half that 
x value, g{2x) oc dF 2 (x, Q 2 )/d(\nQ 2 ) [IB] - The kinematic 
range of the world data for the gluon distribution in nu- 
clear targets is shown in the right panel of Fig.[2]as used in 
a recent global analysis of nuclear modifications to PDFs 
|H1I7]. Such input to the nuclear gluon density is available 
only for x > 0.02 because nuclear DIS has been restricted 
to fixed target experiments. As will be discussed below, 
the study of d(p)+Au collisions at y/s]vjv=200 GeV can 
provide direct sensitivity to the nuclear modification of the 
gluon density for x values on the order of x ~ 0.001 and 
can test the universality of the nuclear gluon density in 
the range 0.02 < x < 0.1. Comparable sensitivity in DIS 
to such low x would require an electron-ion collider |T3] . 
Measuring the PDF with quark and gluon probes allows 
us to get to x and Q 2 values where saturation phenomena 
might be present. 

3 Tests of parton saturation 

Factorization theorems allow us to add cross sections rather 
than quantum amplitudes, with partons considered quan- 
tum mechanically independent of each another. Within 
this picture, we are tempted to imagine that the gluon 
distribution of a nucleus might be obtained by adding the 
gluon distributions for each nucleon, with some account- 
ing for the relative motion of the nucleons in the nucleus. 
While perhaps true for large x processes, at small x the un- 
certainty principle tells us that the partons will all overlap 
in the longitudinal direction, so the partons do not inter- 
act independently. The front surface partons will interfere 
or shadow the back surface partons of the nucleus. For 
more than 20 years it has been recognized |22 that the 
quantum independence of partons cannot extend to very 
small x where the gluon density is very large. 
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Fig. 3. Diagram showing the boundary between possible 
"phase" regions in the r = ln(l/a;) versus ln(Q 2 ) plane |21|. 
BFKL evolution |22| results in expected exponential growth 
of the gluon density at fixed Q 2 and increasing r. DGLAP 
evolution is discussed in the text. 



To determine the scale at which collective behavior 
might become evident, the uncertainty principle suggests 
that a scattering process at fixed pr will probe a trans- 
verse area approximately given by S(pr) — t:{Ti/pt) 2 ■ For 
example at pr = 2 GeV/c, this corresponds to about 0.3 
mb, small in comparison to the proton cross sectional area 
of about 30 mb. The number of gluons that are present and 
that could shadow one another above x is nominally given 
by n g i uons (x) = l x g(x')dx'. At ir=0.01 n gluons « 7 , in- 
creasing by 7-8 for each order of magnitude decrease in the 
lower limit of x. At £=0.001, the product of cross section 
times number of gluons is 5(2 GeV/c) x n g i UO ns (0.001) « 
5 mb. This estimate suggests that for events with these 
kinematics, the chances of finding more than one gluon 
within the transverse resolution of the scattering probe is 
less than 20%. However in a nucleus of mass number A, 
the area of the nucleus grows roughly as A 2 ! 3 while the 
number of gluons would nominally grow proportionally to 
A. Thus, the transverse density of nucleons should grow 
by a factor like A 1 / 3 . For Au, this factor of 6 in transverse 
density suggests that shadowing could be substantial. By 
cc=0.0001, it could become dominant. Of course, at lower 
Pt, the effects would show up at a larger value of x. Real 
predictions for the onset of shadowing vary with the model 
used |23j , but whether shadowing modifies the gluon inter- 
actions at RHIC is an experimental question which must 
be answered with data. 

For p+p collisions at -^5=200 GeV, unlike at lower y/s 
|24j . next-to- leading order (NLO) pQCD calculations |2"5] 
quantitatively describe inclusive particle production down 
to pt of ~2 GeV/c using PDFs and fragmentation 
functions |26U27| that describe the hadronization of the 
scattered partons. Furthermore, di-hadron azimuthal cor- 
relations have the same structure at these moderate pt 
as they do at the highest possible px values, consistent 
with the idea that elastic scattering of quarks and gluons 
is responsible for the particle production |28| . 

Recent measurements at STAR using the prototype 
FPD already indicate that the factorized leading twist 
pQCD calculations work quite well to predict the p + p — ► 
tt°+X cross section in the 3 < r\ < 4 region jTHj. This gives 
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Fig. 4. (Left) pQCD calculation of it — tv° production at large r\ in p+p and d+Au collisions at y / sjvjv=200 GeV |33| . The 
distributions integrate to da/dpr in units of pb/GeV. The smallest x values are probed when tv° — ir° pairs are detected at 
large rj. (Middle) PYTHIA |34| simulation for tv° — tv° production at large r\ in p+p collisions at y/s—2QQ GeV. The rj of the 
associated n° is strongly correlated with the x value of the soft parton involved in the partonic scattering. (Right) HIJING 
simulation for 7r° — n° production at large r\ in d+Au collisions at -^5=200 GeV. Compared to the p+p simulations, the peaks in 
Acj) — — cf> 7T 2 corresponding to elastic parton scattering, sit atop a background from other mechanisms for particle production. 



confidence in the interpretation that at y/s =200 GeV the 
particle production process is dominated by leading twist 
quark- gluon scattering. With the FMS focusing on 7r°— 7r° 
pairs, we will select the low-x component shown by pQCD 
calculations (Fig. to make only small contributions to 
the inclusive measurement. The low-x component of the 
forward pion yield is where shadowing effects are expected 
to be most important [22] • In the middle panel of Fig.0]we 
see that when triggering on a ir° in the range 3 < 77 < 4, 
the rapidity of the second 7r° will reflect the x of the struck 
gluon. The right panel of Fig. 01 shows that elastic par- 
ton scattering is identified above physics backgrounds in 
d+Au collisions. 

There has been considerable recent interest among the 
experts in the application of pQCD in reconciling the 
meaning of shadowing with the idea of universal (factor- 
izable) PDFs. What has emerged recently [3U| is a better 
understanding of just what the universal parton density 
means at small x. The present understanding is that the 
"shadowed" small x distributions should be universal but 
do not strictly reflect the probability for finding a parton. 
Included in the universal factorized functions are built-in 
final-state correlations with other gluons in the proton. We 
now see that even from the strict, leading twist perspec- 
tive, low-x perturbation theory has a different interpreta- 
tion from large x because it always involves a sampling 
of the macroscopic gluon fields. In light of this, there is 
real excitement that a variety of low-x phenomena from 
shadowing to large transverse spin asymmetries may be 
tied together with pQCD in ways not before appreciated. 

Measurements of gluon shadowing at RHIC and the 
LHC will be essential input for models that predict the re- 
lationships between quark distributions and macroscopic 
gluon fields. Among the descriptions of shadowing or sat- 
uration effects is the Color Glass Condensate (CGC) 20, 
I31II32| . an effective-field theory for understanding parton 
saturation. In the CGC picture, the saturation effects are 
associated with a new phase of the gluon field. The onset 
of this phase can be probed by measurements at small x 



and at small Q (related to the produced parton mass and 
the px associated with the scattering). Mapping out the 
boundaries (Fig. [3J for saturation signatures for back-to- 
back jet correlations as a function of x and pt is a primary 
mission of the FMS. 



4 Early results from RHIC 

The first Au+Au collision runs at ^/snn—^OQ GeV re- 
sulted in the observation that high-p-p particle production 
at midrapidity was suppressed relative to expectations re- 
sulting from the scaling of yields from p+p collisions. The 
data also showed that two-particle correlations were sup- 
pressed when the particles were back-to-back (away-side, 
Acp ss 180°) but not when they were fragments of the 
same jet (near-side, Acj) ss or 27r) These observa- 
tions are consistent with a prediction based on radiative 
energy loss of a high-py parton passing through a quark- 
gluon plasma 37 . A d+Au run at y / ijv r /v =200 GeV was 
scheduled early in the RHIC program to eliminate the pos- 
sibility that this was due to initial state effects. For d+Au 
collisions, midrapidity particle production was found to 
have a small enhancement, consistent with the Cronin ef- 
fect and back-to-back correlations [3^1 more closely 
resembled results from p+p collisions than from Au+Au 
collisions. Hence, the suppression of back-to-back correla- 
tions in Au+Au collisions was attributed to the strongly 
interacting matter formed in those collisions, matter that 
was formed by interactions of the low-x gluons. 

Measurements made by the BRAHMS collaboration 
for d+Au collisions [5j showed that inclusive particle pro- 
duction was suppressed as the rapidity of the observed 
particles increased. This provided a hint that the gluon 
distribution in the Au nucleus may be depleted at low-x. 
It is easy to understand how this suppression can occur 
within the standard pQCD picture of particle production. 
In that picture, the quarks and gluons each carry a fraction 
of their parent hadron momentum given by x. They elasti- 
cally scatter and then fragment to the final state hadrons 
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p+p -> 7T° + X, v/s = 200 GeV 




E, (GeV) x„ x F 

Fig. 5. (Left) Invariant cross sections for inclusive ty° production at large rapidity in p+p collisions at -^5=200 GeV |1UI41| 
compared to NLO pQCD calculations |25I2KI27) . (Middle) The data are compared with predictions from PYTHIA H3]. (Right) 
The analyzing power for ir° production at < r\ >=3.8 in p+p collisions at -^5=200 GeV |10|. The curves are predictions from 
pQCD models 1461481501 evaluated at p T =L5 GeV/c. 



observed with a given pt and at a given rj. For collinear 
parton pairs, it is easily shown that 

where the Feynman-x variable is xp = ^Pl/Vs, relevant in 
the limit rji >> 772, and pi is the longitudinal momentum 
component of the large 77 particle. For inclusive particle 
production, one of the two jets, or its hadronic surrogates, 
is observed at 771 and the second jet has a broad 772 distri- 
bution and A(f> w 180°. By detecting a high energy hadron 
at large 771 , initial states with a large- a; parton (most prob- 
ably a quark) and a low-x parton (most probably a gluon) 
are selected. For each unit rapidity increase, the average 
x of the gluon from the initial state parent hadron is de- 
creased by e. A similar decrease in x is a consequence of 
studying particle production in collisions at higher sfs. 
Hence, the observed suppression of particle production at 
increasing rapidity can be interpreted as a reduction in 
the probability of finding gluons in the nucleus at small x. 

The BRAHMS results |S] were confirmed in measure- 
ments by PHENIX gOj. The STAR collaboration also made 
measurements of large rapidity particle production and 
produced a limited data sample for two-particle correla- 
tions involving a large rapidity ir° |28| . measured with 
the FPD. The topology of the energy deposition in the 
FPD allows for robust measurements of the energy and 
direction of neutral pions. Since the 7r° is a pseudoscalar 
particle, kinematic distributions of its diphoton decay are 
exactly calculable in any frame of reference. Hence, cali- 
brations of the FPD response can be obtained at the level 
of ~1% simply by requiring a fully consistent response of 
all cells of the calorimeter to the photons produced by 
the 7t° — > 77 decay. This same technique will be used for 
the FMS. Figure also shows that the simulation code 
PYTHIA [31] is able to reproduce the absolute cross sec- 
tion of the produced 7T° from p+p collisions. In addition, 
nearly all features of the A<j> distributions for charged 



hadrons with \rj\ < 0.75 coincident with forward 7r° are 
reproduced by PYTHIA. This gives us a tool to guide the 
design of the FMS and the interpretation of the data. 



5 Proton spin with the FMS 

Our understanding of the two-particle correlations involv- 
ing large rapidity particles, and our ability to use these 
correlations to measure nuclear gluon distributions, is a 
direct result of the methodology developed to understand 
the first spin asymmetry measurements at RHIC. It is 
no surprise that the FMS will also be a powerful tool in 
studying the spin structure of the proton. 

An early prediction of pQCD was that, at leading twist 
and with collinear factorization, the chiral properties of 
the theory would make the analyzing power (A n) small 
for particles produced with transversely polarized proton 
beams |42| . An is derived from the azimuthal asymme- 
try of particle yields from a transversely polarized beam 
incident on an unpolarized target. However, from AGS en- 
ergies 03] to Fermi Lab energies @3] and most recently at 
STAR JO] (Fig. EJ), a large transverse single spin asym- 
metry has been observed. The consistent trend is that An 
in p j + p — > ir° + X increases rapidly for x f above about 
0.3. Transverse single spin asymmetries have also been 
observed in semi-inclusive DIS from polarized targets 0H] 
and experimental studies of these spin effects is an ac- 
tive area of research. The FMS is ideally suited to extend 
these studies. Calculations of twist-3 contributions 0H] to 
the observed An provide terms that may be related to 
macroscopic gluon fields in the polarized nucleon. 

There are multiple phenomenological effects that have 
been identified as possible sources for the large An, but 
only two that are expected to be large. One is the Sivers 
effect |47U48| , which is an initial state correlation between 
the parton intrinsic transverse momentum kx and the 
transverse spin of the nucleon, An oc St ■ (P x kx), where 
P is the beam momentum and St is the transverse pro- 
ton spin. In the Sivers framework, An is sensitive to the 
contribution of quark orbital angular momentum to the 
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nucleon spin. Large An is the result of a spin dependent 
Pt trigger bias favoring events where kx is in the same 
direction as px- 

If the Sivers effect is present, we can further character- 
ize it with a measurement of the away side jet. The kx of 
the initial-state parton can be measured when final-state 
jet pairs are not exactly back to back (Acf> = 180°) |55| . 
The spin dependence of this kx measurement is exactly 
what the Sivers model predicts. 

While the Sivers effect connects An to the orbital an- 
gular momentum of quarks, a second effect, called the 
Collins effect |49H50| . is directly sensitive to the transver- 
sity distribution function , related to the transverse 

polarization of quarks (and antiquarks) in a transversly 
polarized proton |53| . Here the quark scatters, preserv- 
ing its transverse spin, and then fragments into pions and 
other hadrons. The fragmentation function reveals the po- 
larization of the fragmenting quark and thus the initial 
quark state. In this example, the asymmetry does not ap- 
pear in the jet production directly, but only in the frag- 
mentation. The jet axis would not show the transverse 
asymmetry, but a pion fragment would. Recent calcula- 
tions that include the full kx dependence in the convo- 
lution integrals provide some indication that the Collins 
effect may be small |54|. 

The FMS will be able to distinguish between these 
mechanisms. By looking at pairs of same-side neutral pi- 
ons, we can measure An as a function of the two pion 
kinematics. With the FMS we will separately measure the 
contributions to An that comes from the jet axis versus 
that which comes from the jet structure. Many theory pa- 
pers have studied this problem: however, the need for data 
is great. The FMS STAR experiments on transverse po- 
larization will provide to theorists the necessary input to 
determine the relative contributions from the Sivers effect 
and the Collins effect. 



6 FMS Configuration 

The STAR FPD now taking data acts as the prototype 
for the proposed FMS. Our results with the FPD demon- 
strate the feasibility of large rapidity measurements with 
electromagnetic calorimetry in both p+p and d+Au col- 
lisions at v /s]^v=200 GeV at RHIC. Each FPD calorime- 
ter is a 7 x 7 matrix of 3.8cm x 3.8cm x 45cm lead-glass 
cells that can be positioned in the range 3.3 < 77 < 4.0. 
These are identical to the 684 small cells to be used in 
the FMS and the techniques we have developed for FPD 
tuning and analysis are directly applicable to the FMS. 
By implementing the FMS in STAR, the study of 7 - tt° 
and 7T° — 7T° correlations in both rapidity and azimuthal 
angles is enabled over — 1 < 77 < 4. The FMS also al- 
lows inclusive measurements of ir° , 7 and the production 
of heavy mesons that decay to all 7 final states over a 
broad range of rj and px- A schematic of the proposed 
FMS detector is shown in Fig. H3 along with examples of 
the mass resolution from a topological analysis [SI] of the 
energy deposition for data obtained with the current FPD. 
We expect to achieve 1% accuracy in the calibration of the 



cH-Au->n°+X, Vs„=200 GeV, <77> = 4.0 
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Fig. 6. (Top) Di-photon invariant mass distributions from 
the STAR FPD for d+Au collisions at ^s^iv=200 GeV at 
< 77 >=4.0 relative to the incident deuteron beam. (Bottom) 
Layout of the proposed FMS. 



FMS, just as we have done for the FPD. Energy resolution 
of < has been demonstrated. Simulation studies 

compare well to FPD data and show that we can expect to 
locate the tt° to better than 0.5 cm (RMS). Our expected 
7T° mass resolution is ^23 MeV/c 2 based on experience 
with the FPD. For neutral pions with 20 < E < 60 GeV 
the reconstruction efficiency of the FPD is just given by 
geometry. To predict event yields, we used a very conser- 
vative estimate of 35% for reconstruction and geometric 
efficiency The FMS will provide complete azimuthal cov- 
erage for the pseudorapidity interval 2.5 < 77 < 4.0 and 
will be built from existing lead-glass cells. 



7 Plan for measurements 

We will concentrate on measurements directed at our three 
immediate physics goals for d+Au and for p+p running. 
In a future d+Au run at y / s=200 GeV we will determine 
the gluon distribution in the gold nucleus. The large ac- 
ceptance of the FMS will provide good statistics for com- 
parisons of identified ir° yields in p+p and d+Au colli- 
sions over a very broad range in px and pseudorapidity. 
We will measure the correlations between a trigger parti- 
cle in the FMS (either a 7r° or a 7) and a second particle 
(jet surrogate) in the TPC, BEMC, EEMC or FMS. The 
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threshold pt for the trigger and coincident particles will 
be investigated over the range from 1 < px < 4 GeV/c. 
The large size and high granularity of the FMS allow us 
to identify neutral pions at energies down to ~10 GeV, 
where hadronic energy deposition becomes a significant 
background. At maximum rapidity, this determines our 
lowest observable px threshold. As seen in Fig. the r\ of 
the coincident particle correlates with the gluon x value. 
The Atfi distribution reveals the scattering history. To de- 
velop the full picture of the gluon distributions, we will 
investigate the dependence of tt° — ir° correlations on Q 2 . 
To determine the dynamical origin (Sivers or Collins) of 
the observed transverse single spin asymmetry for forward 
7T° production, we will make similar correlation measure- 
ments in the polarized p+p run. 

As shown in Fig. pQCD calculations [33] suggest 
that the lowest x values for the gluon density are probed 
when both jets from the elastic parton scattering are pro- 
duced at large rapidity. Both ir° + ir° and ir° + 7 final 
states can be analyzed in the FMS. The pr of the ir° 
must be large enough to favor elastic parton scattering 
over inelastic scattering [2Hj although both contributions 
are contained in PYTHIA and NLO pQCD calculations. 
The large size of the proposed FMS enables the use of 
isolation cuts to distinguish between ir° decay photons 
and direct photons. The dominant subprocess for direct 
photon production at RHIC is QCD Compton scattering, 
qg — > jq. This can also be used to probe the small x gluon 
density with only minimal physics backgrounds. 

FigureHshows a HIJING 1.381 [35] simulation of 7.5 x 
10 s minimum-bias d+Au events at ^/snn— 200 GeV, from 
which all ir° — 7r° pairs with the specified pt and r\ are se- 
lected and used to compute A(j>. Unlike the case for p+p 
collisions, the elastic parton scattering peaks in the A<p 
distribution sit atop a background from the nuclear colli- 
sion. Despite the background, the elastic parton scattering 
is readily discriminated and can be identified as the ex- 
pected peak in A<j>. We can quantitatively describe these 
distributions from 1 < A<f> < 5.28 radians by a Gaus- 
sian function, used to model the peak, and a constant 
background. A best fit to the A<j) distribution in Fig. [4] 
results in 2.3 x 10 3 events in the peak for simulations 
done without shadowing and ~ 1.8 x 10 3 events in the 
peak for the simulations done with shadowing. We can 
expect that a 10-week d+Au run will allow us to sam- 
ple > 6 x 10 10 minimum-bias interactions, based on RHIC 
performance for d+Au collisions achieved in the last weeks 
of run 3. Accounting for detector efficiencies for the FMS 
and endcap EMC, the simulations suggest we will observe 
at least 8 x 10 3 events in the A<f> peak without shadowing. 
This is enough to investigate the spatial dependence of 
the nuclear gluon density [T^ITH] using particle multiplic- 
ity measurements in the Au beam direction made by other 
STAR subsystems to determine sensitivity to the impact 
parameter of the collision. Broadening or disappearance 
of the away-side correlation could signal the transition to 
a macroscopic gluon field at sufficiently low x. 

For the polarized proton running, we base estimates 
of the rate of near-side ir° — ir° pairs on PYTHIA |34) . 



The STAR measurements (Fig. show that Am is small 
below if ~ 0.4, and increases monotonically as xf of 
the forward tt° increases. In the Sivers picture, Am should 
be associated with the forward jet and should be present 
for 7T° — 7T° pairs from the same jet. We would expect a 
large Am when xpx + %F2 > OA. When one tt° is ob- 
served at 3 < rj < 4 with xpi > 0.25 and a second ir° 
is observed with ij < 4 and Xfi > 0.15 and the 7r° — ir° 
pair has \rji — 772 1 < 0.5, the simulated A<fi correlation 
shows a jet-like near-side correlation peak sitting atop a 
uniform background attributed to the underlying event. 
We expect ~ 1.5 x 10 4 tt° — tt° events in the near-side 
jet-like peak for 1 pb _1 of integrated luminosity for po- 
larized p+p collisions at y/s = 200 GeV with the FMS. 
For a beam polarization of 50% this would result in a 
statistical error on the analyzing power of 5 Am ~ 0.01. 
The Collins mechanism attributes Am to the correlation 
between the momenta of two hadrons from the same jet 
and the proton spin vector. The transverse momentum 
associated with jet fragmentation producing a 7r° with 
xfi > OA and 3 < 77 < 4 can be measured by detect- 
ing a second ir° with 77 < 4 and Xfi > 0.15 and requiring 
the 7T° — 7T° pair have |ryi — 772 1 < 0.5. Again, the A<j) cor- 
relation shows a jet-like near-side correlation peak sitting 
atop a uniform background. For these kinematics, we ex- 
pect 2 x 10 3 7T° — 7T° pairs in the jet-like near-side A<p peak 
for lpb -1 of integrated luminosity. If a non-zero Collins 
effect is observed, then larger integrated luminosity sam- 
ples would be required to map out the x dependence of 
the transversity structure function. 

8 Summary 

Early experimental results from RHIC, coupled with the- 
oretical developments that demonstrate the robustness of 
NLO pQCD calculations at collider energies down to mod- 
erate pt values, suggest that a quantitative determination 
of the gluon density in a heavy nucleus can be obtained 
at an order of magnitude lower x than available from DIS 
on nuclear targets. Such studies require improved instru- 
mentation in the forward direction at RHIC. In addition 
to providing crucial information to understand the initial 
state of heavy-ion collisions that may lead to a quark- 
gluon plasma, future low- a; studies at RHIC may establish 
the existence of a macroscopic gluon field. Improved for- 
ward instrumentation can also disentangle the dynamical 
origin of transverse single spin asymmetries. 
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